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ABSTRACT 

The physical modeling of active regions (ARs) and of the global coronal is receiving increasing 
interest lately. Recent attempts to model ARs using static equilibrium models were quite successful 
in reproducing AR images of hot soft X-ray (SXR) loops. They however failed to predict the bright 
EUV warm loops permeating ARs: the synthetic images were dominated by intense footpoint emission. 
We demonstrate that this failure is due to the very weak dependence of loop temperature on loop 
length which cannot simultaneously account for both hot and warm loops in the same AR. We then 
consider time-dependent AR models based on nanoflare heating. We demonstrate that such models 
can simultaneously reproduce EUV and SXR loops in ARs. Moreover, they predict radial intensity 
variations consistent with the localized core and extended emissions in SXR and EUV AR observations 
respectively. We finally show how the AR morphology can be used as a gauge of the properties 
(duration, energy, spatial dependence, repetition time) of the impulsive heating. 

Subject headings: hydrodynamics; Sun : Coronal Heating, Sun : Corona 


1. INTRODUCTION 

Modeling ensembles of coronal loops in active regions 
or over the full Sun is a rapidly emerging new field in 
the study of the solar corona. Such studies attempt to 
reproduce generic properties of ARs and of the global 
corona such as space integrated intensities and overall 
morphology. This allows for the determination various 
properties of coronal heating such as its dependence on 
loop length and magnetic field thereby testing coronal 
heating mechanisms. Moreover such studies eliminate 
potential selection effects which may enter into studies 
of individual loops. Finally, AR and global coronal mod- 
els pave the way for the construction of physics-based 
models of the EUV and SXR irradiance, an important 
contributor to Space Weather conditions. 

The majority of the attempted AR and global coro- 
nal models were based on static equilibrium coronal loop 
models produced with steady heating (Schrijver et al. 
2004; Mok et al. 2005; Warren & Wiriebarger 2006; 
Lundquist, Fisher & McTiernan 2007). They were quite 
successful at reproducing the general appearance of the 
corona in hot emissions (> 3 MK) observed in the SXR 
by the Soft X-ray Telescope (SXT) on Yohkoh. However, 
synthetic EUV images from the static models did not 
show evidence of the warm loops seen in the real obser- 
vations, but instead were dominated by intense footpoint 
(moss) emission. Steady heating that is highly concen- 
trated near the coronal base can produce time dependent 
behavior called thermal non-equilibrium. This gives rise 
to transient EUV loops (Mok et al. 2008; Klimchuk, 
Karpen, & Patsourakos 2007), but whether their proper- 
ties are fully consistent with observations has yet to be 
determined. 

Most theories of coronal heating predict that energy 
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is released impulsively on individual flux strands (Klim- 
chuk 2006). This includes both AC (wave) and DC 
(reconnection-type) heating. Following this paradigm 
coronal loops are viewed as ensembles of unresolved, im- 
pulsively heated strands. A very promising idea first pro- 
posed by Parker (1988) is that the coronal field becomes 
tangled on small scales due to the random footpoint mo- 
tions associated with turbulent photospheric convection. 
Current sheets develop at the interfaces between individ- 
ual misaligned strands, and it has recently been shown 
that an explosive instability called the secondary insta- 
bility occurs when the misalignment angle reaches a crit- 
ical value (Dahlburg, Klimchuk & Antiochos 2005). This 
may be the physical nature of the nanoflares postulated 
by Parker. The heating function we have used in our sim- 
ulations is appropriate to nanoflares that are initiated at 
a critical angle. We note that the upward Poynting flux 
associated with observed photospheric field strengths and 
observed photospheric velocities is consistent with coro- 
nal heating requirements (e.g. Abramenko, Pevtsov & 
Romano 2006), but the efficiency of field line tangling is 
only now being addressed quantitatively with high reso- 
lution magnetogram movies from Hinode. We also note 
that direct observations of the kinds of nanoflares we 
are discussing are scarce (Katsukawa & Tsuneta 2001). 
The small distinct brightenings that are sometimes called 
nanoflares are much different from the unresolved energy 
releases on long field lines that we consider here. 

But why do static equilibrium models fail to reproduce 
the coronal emission patterns in the EUV? Can impul- 
sive heating models do any better? And how can AR 
morphology be used to infer the properties of impulsive 
heating? With this paper we address these important 
questions. 

2. AR SIMULATIONS 

For the calculations reported in this paper we use the 
0D hydrodynamic loop model we call enthalpy-based 
thermal evolution of loops (EBTEL), described in Klim- 


chuk, Patsourakos k Cargill (2008). For a given tempo- 
ral profile of the heating EBTEL calculates the temporal 
evolution of the spatially averaged temperature, density, 
and pressure along the loop. It also provides the Dif- 
ferential Emission Measure distribution (DEM(T)) for 
both the coronal and footpoint (i.e., transition region) 
sections of the loop. For steady heating, static equilib- 
rium solutions are calculated. EBTEL is capable of rel- 
atively accurately mimicking complex ID hydrodynamic 
simulations, with however much less demands in terms 
of CPU time, which makes it particularly useful for para- 
metric investigations of multitudes of loops. We calcu- 
lated hydrodynamic models for 26 loops with lengths in 
the range 50-150 Mm, typical of observed AR loops. For 
the construction of the AR images we assumed that the 
loops have semi-circular shapes and are nested the one 
on top of the other, forming an arcade which emulates 
the simplest form of AR topology. 

2.1. Static Models 

We first calculated a static equilibrium AR model. The 
steady volumetric heating H supplied to a loop with 
length L was : 

H = H 0 (L/L 0 ) a , (1) 

with Ho the heating magnitude, Lq the length of the 
shortest loop and a a scaling-law index which depends 
on the details of the specific coronal heating mechanism 
(e.g., Mandrini et a!. 2000). We chose Lq = 50 Mm, 
// o —0.01 erg cirrus V and a = —2.8. This particu- 
lar a corresponds to heating associated with the tan- 
gling of the magnetic field by photospheric convection. 
A nanoflare occurs when the misalignment between ad- 
jacent flux strands reaches a critical angle. Similar a val- 
ues were found to provide the best match between static 
models and AR and full Sun SXR images (Schrijver et al. 
2002; Warren k Winebarger 2006; Lundquist, Fisher k 
McTiernan 2008). In concert with these studies, H 0 was 
chosen so that the temperature of the shortest and con- 
sequently hottest (according to Equation 1) loop of the 
arcade was « 4 MK, consistent with the temperature of 
bright SXR loops in AR cores. The DEM distributions 
for the coronal and footpoint sections of each simulated 
loop were folded with the temperature response functions 
of the 171 A and AlMg channels of TRACE and SXT 
respectively to determine the corresponding intensities. 
For the calculation of the intensities we assumed that 
simulated loops had a diameter of 3 Mm, consistent with 
typical widths of observed loops. 

In Figure 1 (left column) we plot the variation of in- 
tensity across the loop arcade while in Figure 2 (left col- 
umn) we display the corresponding synthetic images. For 
building the images we assume that the AR is viewed 
face-on and for clarity reasons we display every fourth 
loop of our arcade. We applied convenient box-cars to 
the images to emulate the different spatial resolution of 
SXT (« 5 arcsec) and TRACE (« 1 arcsec). We assumed 
that the coronal emission is distributed uniformly along 
the loop. This is reasonable because our loops are shorter 
than a gravitational scale height, so there will be mini- 
mal gravitational stratification. Furthermore, the coro- 
nal temperature varies by only about 50% along most of 
the length of equilibrium loops (Klimchuk et al. 2008). 
Temperature and density of course vary dramatically in 


the transition region, but the thickness of transition re- 
gion is generally unresolvable, so we spread the emission 
over 2 Mm for convenience and clarity. The magnitude 
of the integrated transition region emission is correct. 

From Figure 1 we note that the SXT emission is con- 
siderably weaker at the footpoints than in the corona. 
This is because the footpoint temperatures are generally 
below 2 MK, where SXT has greatly reduced sensitivity. 
Furthermore, from Figure 2 the synthetic TRACE image 
is completely dominated by the footoints. There is little 
evidence of EUV loops, which is at odds with the mul- 
titudes of EUV loops seen in the majority of observed 
ARs. Our results are broadly consistent with previous 
studies employing static heating. 

Why is it that static equilibrium models fail to predict 
bright EUV loops together with bright SXR loops? This 
’’pathology” is related to the fact that under static equi- 
librium conditions the loop temperature has a very weak 
dependence on the loop length. For instance, the Rosner 
et al. (1978) scaling law predicts that the apex temper- 
ature is related to the heating rate and length according 
to T a oc // 2/ ' 7 /// 7 . Therefore, for the employed value 
a = —2.8 we have from Equation 1 that T a oc L~ 0 - 22 . 
This means that the temperature is reduced by a factor of 
only 3 0 22 = 0.8 in going from the shortest loop (4 MK) 
to the longest loop (3.1 MK) in the arcade. As shown 
Figure 3, our simulations closely follow the above scaling 
law. Therefore, none of the loops has 1-2 MK plasma in 
the coronal section, which is necessary to produce sig- 
nificant EUV emission. All of this plasma resides at the 
footpoints, which is where the strong TRACE emission 
originates. Obviously, we could decrease Hq to produce 
strong TRACE emission in the corona, but then the SXT 
emission would be dramatically reduced. The only way 
to have both bright TRACE loops and bright SXT loops 
in the same arcade is for the heating rate to have a much 
stronger dependence on loop length. Our example ar- 
cade requires a ~ -6. We are not aware of any coronal 
heating mechanism with such an extreme dependence on 
L (e.g. Mandrini et al. (2001)). 

2.2. Impulsive Models 

We then considered models with impulsive heating, us- 
ing the same loop arcade of the previous section. We 
started with static equilibria having an average coronal 
temperature near 0.5 MK. We heated the loops with a 
triangular pulse lasting 50 s and let the loops cool for 
8500 s, by which time the temperature had cooled be- 
low 1 MK. The amplitude of the heat pulses varied from 
loop to loop according to Equation 1, with a = —2.8, 
as before, and Ho = 2.5 erg cm" 3 s -1 . This produced an 
average DEM for the arcade that peaks near 2.0 MK, 
similar to that of observed active regions (e.g., Brosius 
et al. 1996). Each loop reaches a maximum tempera- 
ture exceeding 30 MK (Figure 3), but this happens very 
early in the heating event, when the density is very low. 
The DEM- weighted mean temperature Tdbm is only 
1-2 MK. 

We determined temporally averaged TRACE and SXT 
intensities for each loop simulation. Time averaging over 
the duration of a loop simulation is equivalent to taking 
a snapshot of a loop containing a large number of im- 
pulsively heated unresolved strands at different stages of 
heating and cooling. As before, we produced synthetic 


images of the arcade by assuming that the time-averaged 
intensities are uniform along each the loop. Our ID hy- 
dro simulations indicate that this is a reasonable approx- 
imation (e.g., Klimchuk et al. 2006). The right columns 
of Figures 1 and 2 show the intensity variation across the 
arcade and the synthetic image, respectively. 

We note in Figure 1 that the TRACE emission from 
the footpoints is a factor of « 3-100 smaller for impulsive 
heating than for static heating, wheres the coronal emis- 
sion is about the same in the two cases. The brightness 
contrast between the footpoints and corona is therefore 
significantly reduced for impulsive heating. This leads 
to a TRACE image in which both the coronal and foot- 
point emissions can be readily discerned (Figure 2 right 
panel). This is not true for the static model (left panel), 
where the coronal emission is overwhelmed by the foot- 
point emission. The footpoint to corona intensity ra- 
tios I foot/ 1 corona are of order 10 in the impulsive model 
and 1000 in the static model. Observed values are in 
the range « 2-20. The smaller observed ratios could be 
due to spicular absorption of the footpoint emission (e.g., 
Daw, Deluca & Golub 1995; De Pontieu et al. 1999). Us- 
ing the analytical expressions of Anzer & Heinzel (2005) 
for absorption at TRACE wavelengths and typical phys- 
ical parameters of spicules given in Table- 1 of Tsiropoula 
& Schmieder (1993) we found that attenuation factors of 
the 171 TRACE emission of about 10 can be achieved. 
As a matter of fact a recent SUMER/EIS study of moss 
intensities formed above and below the head of the hy- 
drogen Lyman continuum at 912 A demonstrated that 
indeed sizeble absorption occurs over moss regions (De 
Pontieu et al. 2008). Future AR models would proba- 
bly need to incorporate absorption effects. Warren and 
Winebarger (2007) modeled an observed AR using ID 
simulations and also found that impulsive heating in- 
creases the visibility of EUV loops compared to static 
equilibrium. 

It can also be seen in Figure 1 that brightness of the 
corona in TRACE relative to SXT is roughly two orders 
of magnitude larger in the impulsive model than in the 
static model. This suggest that bright TRACE loops are 
much more likely to be seen together with bright SXT 
loops in the same active region if the heating is impulsive. 

The increased TRACE-to-SXT coronal intensity ratio 
and the reduced footpoint-to-corona TRACE contrast 
can both be understood as follows. In static equilibrium, 
the optimum temperature for a particular waveband oc- 
curs either in the corona or at the footpoints, but not 
at both locations. For apex temperatures > 2 MK, the 
coronal plasma is too hot to emit appreciable TRACE 
emission, and only the footpoints are bright. This is 
the case for all of the loops in the static arcade. Impul- 
sively heated loops are fundamentally different in that 
they experience a wide range of coronal and footpoint 
temperatures over the course of their evolution. When 
an impulsively heated loop starts to cool, TRACE emis- 
sion occurs first at the footpoints and then over the full 
length of the loop when the coronal temperature drops 
below 2 MK. A bundle of impulsively heated strands will 
therefore have significant TRACE emission both in the 
corona and at the footpoints. 

A final interesting property revealed in Figure 1 is the 
distribution of coronal emission across the arcade. SXT 
intensities decrease rapidly with loop length, especially in 


the impulsive model, while TRACE intensities decrease 
much more slowly. As a consequence, the SXT emission 
is concentrated in the core of the arcade, and the TRACE 
emission is much more extended (Figure 2). This agrees 
well with AR observations. The strong decrease in the 
SXR emission with increasing L is because Tdem, the 
dominant temperature of the plasma, drops below 2 MK 
for the longer loops (Figure 3). The sensitivity of SXT 
is a rapidly decreasing function of temperature in this 
range. On the other hand, the Tdbm values are in the 
range where TRACE has good sensitivity, so the TRACE 
intensity gradient across the arcade is more shallow. 

We then examined how the AR morphology depends 
on the properties of the impulsive heating. We first con- 
sidered the effect of varying the index a of Equation 1 by 
taking the following values: -4, -3, -2, -1. All other as- 
pects of the simulation were the same as for the a = —2.8 
base simulation discussed above. The resulting intensity 
variations across the loop arcade for TRACE and SXT 
are given in Figure 4. Not surprisingly, we found that 
a stronger dependence of impulsive heating on L (i.e. a 
more negative a) produces a steeper intensity drop-off 
across the arcade. Note that that if we were to observe 
our arcade off-limb, the model with a=-l would yield 
an intensity scale-height consistent with off-limb observa- 
tions of the unresolved EUV corona (Cirtain et al. 2005). 

The impact of the time-interval between successive 
nanoflares, Ar nano , assuming the dependence Ar„ ano = 
Aro(L/Lo) 7 was then examined. We consider two cases: 
(1) 7 = -2 and Ar 0 = 8500 s; (2) 7 = 2 and Ar 0 = 250 
s. For each loop we considered 3 consecutive nanoflares 
separated by Ar„ ano . The nanoflare heating rate and 
duration were the same as in the base simulation which 
implies the same averaged energy per nanoflare. When 
AT nano decreases with L (case 1) the resulting images 
(left column of Figure 5) are characterized by TRACE 
emission concentrated around the AR core while the SXT 
emission is extended. This is because the longest loops in 
the arcade have A T nano of about 850 s which is not suf- 
ficient to let them cool down to TRACE temperatures. 
On the other hand, when A T nano increases with L (case 
2), the resulting AR images are characterized by bright 
core emission in SXT and lack of TRACE loops (right 
column of Figure 5). This is because shorter loops do 
not have the time to cool down to TRACE temperatures 
and they remain in quasi-steady state conditions in SXT 
temperatures, as also shown in Section 2.1. This model 
can explain observations of an AR with SXT loops but 
no TRACE loops ( Antiochos et al. 2003) . 

We then investigated the effect of the magnitude of 
the nanoflare energy on the SXT and TRACE coronal 
and moss intensities respectively. These intensities are 
particularly useful when studying the cores of active re- 
gions and can constrain the properties of the heating 
(e.g. Winebarger, Warren & Falconer 2007). On top 
of the base model described in the first paragraph of 
this section, we considered 3 additional AR models em- 
ploying nanoflare heating magnitudes 2,6 and 10 times 
the one used in the base simulation (see also Equation 
1). All other aspects of these simulations were the same 
as for the base simulation. The variation of the tem- 
porally averaged intensities across the loop arcade for 
the 4 models is shown in Figure 6. What can be seen 
in this Figure is that increasing nanoflare energy leads 


to higher SXT coronal to TRACE moss intensity ra- 
tio Isxr/Ieuv — moss - For instance, for the inner loops 
in the modeled AR (i.e. its core) I S xr/ I suv-moss in- 
creases from ps 1.3 to 4.4 with a 10-fold increase in the 
nanoflare energy. This means that IsxR/lEUV~-mo SS 
tracks in a rather sensitive way the nanoflare energy and 
can be used as its diagnostic, is it worth mentioning that 
the determined trend goes in the direction of decreasing 
the discrepancy with observations that can cause prob- 
lems with static models. Possibly the spicular absorp- 
tion discussed in the previous paragraphs could bring 
Isxr/Ieuv — moss even to a better agreement with the 
observations 

We finally found that the nanoflare duration does not 
alter the AR morphology for both SXT and TRACE. 
This comes to no surprise given both instruments are 
mostly sensitive to the late cooling of the impulsively 
heated loops, when any differences in the loop response 
for different nanoflare durations, had been long ago 
smeared out (e.g. Patsourakos & Klimchuk 2005). 

3. CONCLUSION 

With this work we considered static and impulsive 
heating models of active region arcades. In concert with 
previous investigations, we found that static models can- 
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not simultaneously reproduce bright EUV and SXR loops 
withing the same AR. We showed that this is due to the 
shallow dependence of loop temperature on loop length 
implied by all reasonable coronal heating scenarios. We 
found that impulsive heating models agree much better 
with observations than do static models. In particular, 
impulsive heating produces (1) a reduced brightness con- 
trast between the corona and footpoints in TRACE ob- 
servations, (2) an increased TRACE-to-SXT coronal in- 
tensity ratio, and (3) enhanced SXT emission in the core 
of the active region and extended TRACE emission. We 
finally showed that the AR morphology depends rather 
sensitively on the properties of impulsive heating, like its 
spatial dependence and the time interval between succes- 
sive nanoflares, and can therefore be used to determine 
the properties of the heating. Our study paves the way 
for detailed comparisons between multi-temperature ob- 
servations of ARs and models based on impulsive heating 
and detailed reconstructions of the coronal magnetic field 
from extrapolations. 
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Fig. 1. Variation of the intensity across the loop arcade for static equilibrium (left column) and for impulsive heating (right column) 
for SXT and TRACE. We plot the intensities for both the coronal (solid line) and transition region section (dashed line) of the loops. For 
impulsive heating the time-averaged intensities are plotted. Intensities are in units of DN/pix/s. 
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Fig. 2. — SXT and TRACE images of the simulated AR loop arcade for static equilibrium (left column) and impulsive heating (right 
column). The images for the impulsive heating correspond to a time-average over the corresponding simulations. It is assumed that the 
AR is observed face-on. The images are smoothed with boxcars consistent with the instrument resolution. The AR baseline has a length 
of PS 100 Mm. Intensity represented with color increasing from black to red to white. Each image is normalized individually. 
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Fig. 3. — Variation of the temperature across the loop arcade. Static equilibrium model (solid line) and fit tx L -0 - 22 (dash-triple dot). 
Impulsively heated model: maximum temperature divided by 10 (dashed line); DEM - weighted mean temperature (dashed-dotted line). 
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Fig. 4. Variation of the normalized coronal intensity across the loop arcade for different values of the power-law index a determining 
how the nanoflare energy is distributed across loops with different L, for TRACE (left panel) and SXT (right panel). 
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Fig. 5.— Investigating the effect of the time-interval A r nano between successive nanoflares. Synthetic TRACE and SXT images for 
Arna-no which decreases (increases) with loop length left (right) column. The AR baseline has a length of as 100 Mm. Intensity represented 
with color increasing from black to red to white. Each image is normalized individually. 
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Fig. 6. Variation of the coronal and footpoint intensity across the loop arcade for SXT (left panel) and TRACE (right panel) respectively 
for impulsively heated AR models with different nanoflare energies. The nanoflare magnitude for each AR model are 4 (dots), 6 (long 
dashes), 10 (dash-dot) times that of the base model (solid). The time-averaged intensities are plotted. Intensities are in units of DN /pix/s. 




